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Abstract. We study the effect of non-radial motions on the mass func- 
tion, the VDF and on the shape of clusters of galaxies using the model 
introduced in Del Popolo & Gambera ( 1998a, b; 1999). The mass function 

■ of clusters, obtained using the quoted model, is compared with the sta- 
Q\ \ tistical data by Bahcall & Cen (1992) and Girardi et al. (1998), while the 

VDF is compared with the CfA data by Zabludoff et al. (1993) for local 
clusters. In both cases the model predictions are in good agreement with 

■ the observational data showing once more how non-radial motions can 
reduce many of the discrepancies between CDM model predictions and 
observational data. Besides we study the effect of non-radial motions on 

r-j ! the intrinsic shape of clusters of galaxies showing that non-radial motions 

CLi 1 produce clusters less elongated with respect to CDM model. 
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1. Introduction 

At its appearance the CDM model contributed to obtain a better understanding 
of the origin and evolution of the large scale structure in the Universe (White 
et al. 1987). The principal assumptions of the standard CDM (SCDM) model 
(Liddle & Lyth 1993) are: 

• a flat Universe dominated by weakly interacting elementary particles hav- 
ing low velocity dispersion at early times. The barionic content is deter- 
mined by the standard big bang nucleosynthesis model; 

• critical matter density; 

• expansion rate given by h = 0.5; 

• a scale invariant and adiabatic spectrum with a spectral index, n = 1; 

• the condition required by observations, that the fluctuations in galaxy 
distribution, (Sp/p) g , are larger than the fluctuations in the mass distri- 
bution, (8p/p)p by a factor b > 1. 
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If this last assumption is not introduced, the pairwise velocity dispersion is 
larger then that deduced from observations and the galaxy correlation function 
is steeper than that observed (Davis et al. 1985). After the great success of 
the model in the 80's, a closer inspection of the model has shown a series of 
deficiencies, namely: 

• the strong clustering of rich clusters of galaxies, £ C c( r ) — (r/25ft. _1 Mpc) -2 , 
far in excess of CDM predictions (Bahcall & Soneira 1983); 

• the overproduction of clusters abundance. Clusters abundance is a useful 
test for models of galaxy formation. This is connected to three relevant 
parameters: the mass function, the VDF and the temperature function. 
Using N-body simulations, Jing et al. (1994) studied the mass function 
of rich clusters at z = for the CDM model. They found that the CDM 
model with the COBE normalization produces a temperature function of 
clusters higher than that given by the observations by Edge et al. (1990) 
and by Henry h, Arnaud (1991); 

• the conflict between the normalization of the spectrum of the perturbation 
which is required by different types of observations; 

• the incorrect scale dependence of the galaxy correlation function, £(r), on 
scales 10 100 /i _1 Mpc, having £(r) too little power on the large scales 
compared to the power on smaller scales (Maddox et al. 1990; Saunders 
et al. 1991; Lahav et al. 1989; Peacock 1991; Peacock & Nicholson 1991). 

Several alternative models have been proposed in order to solve the quoted prob- 
lems (Peebles 1984; Valdarnini & Bonometto 1985; Bond et al. 1988; Holtzman 
1989; Efstathiou et al. 1990; Turner 1991; Schaefer 1991; Schaefer & Shafi 1993; 
Holtzman Sz Primack 1993; Bower et al. 1993). Most of them propose in some 
way a modification of the primeval spectrum of perturbations. In two previous 
papers (Del Popolo Sz Gambera 1998a; 1999) we showed how, starting from a 
CDM spectrum and taking into account non-radial motions, at least the problem 
of the clustering of clusters of galaxies and the problem of the X-ray temperature 
can be considerably reduced. 

Here, we study the effect of non-radial motions on the mass function, the VDF 
and on the shapes of galaxy clusters. 

In Sect. 2 we shall use the same model introduced by Del Popolo Sz Gambera 
1998a, b, 1999) to compare the mass function calculated using the CDM model, 
taking into account non-radial motions, with the observed mass function ob- 
tained by Bahcall Sz Cen (1992). Then, we repeat the calculation for the VDF 
and compare the theoretical VDF with the CfA data by Zabludoff et al. (1993). 
In Sect. 4 we study the effect of non-radial motions on the ellipticity of clusters 
and finally in Sect. 5 we give our conclusions. 

2. The mass function and the velocity dispersion function 

One of the most important constraints that a model for large-scale structure 
must overcome is that of predicting the correct number density of clusters of 
galaxies. This constraint is crucial for several reasons (Peebles 1993). The 
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abundance of clusters of galaxies, together with the mass distributions in galaxy 
halos and in rich clusters of galaxies, the peculiar motions of galaxies, the spatial 
structure of the microwave background radiation is one of the most readily 
accessible observables which probes the mass distribution directly. The most 
accurate way of assessing the cluster abundance is via numerical simulations. 
However, there is an excellent analytic alternative, Press & Schechter's theory 
(Press & Schechter 1974; Bond et al. 1991). Press-Schechter's theory states 
that the fraction of mass in gravitationally bound systems larger than a mass, 
M, is given by the fraction of space in which the linearly evolved density field, 
smoothed on the mass scale M, exceeds a threshold S c : 

F(> M) = -erfc ( —=-^ J (1) 

where Rf is the comoving linear scale associated with M. Press-Schechter's 
result predicts that only half of the mass of the Universe ends up in virialized 
objects but in particular cases this problem can be solved (Peacock & Heavens 
1990; Cole 1991; Blanchard et al. 1992). 

The mass variance present in Eq. ([j]) can be obtained once a spectrum, P(k), is 
fixed: 

1 r°° 

a 2 (M) = — , / dkk 2 P(k)W 2 (kR) (2) 
2tt z Jo 

where W{kR) is a top-hat smoothing function: 

W(kR) = -^-^ (smkR- kRcoskR) (3) 

and the power spectrum P(k) oc k n T 2 (k) is fixed giving the transfer function 
T(/c)(here, we adopt that by Klypin et al. 1993), where k is the wave-number 
measured in units of Mpc -1 . This spectrum is valid for k < 30Mpc _1 and 
z < 25. The accuracy of the spectrum is 5%. It is more accurate than Holtz- 
man's (1989) spectrum, used by Jing et al. (1994) and Bartlett & Silk (1993) 
to calculate the mass function and the X-ray temperature function of clusters, 
respectively. The spectrum is lower by 20% on cluster mass scales than Holtz- 
man's (1989). The spectrum was normalized to the COBE quadrupole Q2 = 17yU 
K, corresponding to = 0.66. As shown by Bartlett & Silk (1993) the X-ray 
distribution function, obtained using a standard CDM spectrum, over-produces 
the clusters abundances data obtained by Henry & Arnaud (1991) and by Edge 
et al. (1990). This has lead some authors (White et al. 1993) to cite the clus- 
ter abundance as one of the strongest pieces of evidence against the standard 
CDM model when the model is normalized so as to reproduce the microwave 
background anisotropies as seen by the COBE satellite (Bennett et al. 1996). 
The discrepancy can be reduced taking into account the non-radial motions that 
originate when a cluster reaches the non-linear regime as follows. A fundamental 
role in Press-Schechter's theory is played by the value of 5 C . Using a top-hat 
window function 5 C = 1.7 ± 0.1 while for a Gaussian window the threshold is 
significantly lower. In a non-spherical context the situation is more compli- 
cated. Considering the collapse along all the three axes the threshold is higher, 
whereas the collapse along the first axis (pancake formation) or the first two 
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axes (filament formation) corresponds to a lower threshold (Monaco 1995). The 
threshold, <5 C , does not depend on the background cosmology. 
As shown by Del Popolo & Gambera (1998a; 1999), if non-radial motions are 
taken into account, the threshold S c is not constant but is function of mass, M 
(Del Popolo & Gambera 1998a; 1999): 



1 + 



rta r ta L 2 • dr 
, GM 3 r 3 



(4) 



where <5 co = 1.68 is the critical threshold for a spherical model, r\ is the initial 
radius, r ta is the turn-around radius and L the angular momentum. In terms 
of the Hubble constant, Hq, the density parameter at current epoch, Qq, the 
expansion parameter a and the mean fractional density excess inside a shell of a 
given radius, 5 Eq. (||) can be written as (Del Popolo & Gambera 1998a; 1999): 



i + 



8G 2 
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(5) 



where a is the expansion parameter, and a m ; n its value corresponding to r\. The 
mass dependence of the threshold parameter, 5 c {y), was obtained in the same 
way sketched in Del Popolo & Gambera (1999). The result of the calculation is 




Figure 1. Cumulative mass function calculated using a CDM model 
without taking into account non-radial motions (dashed line) and tak- 
ing account non-radial motions (solid line) compared with Bahcall & 
Cen (1992) data (full dots) and with that of Girardi et al. (1998) (open 
squares) . 



shown in Fig. 1. Here, the mass function of clusters, derived using a CDM model 
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with f^o = 1) h = 1/2 normalized to Qcobe = 17/t K and taking into account 
non-radial motions (solid line), is compared with the statistical data by Bahcall 
& Cen (1992) (full dots) and with that of Girardi et al. (1998) (open squares) 
and with a pure CDM model with Oo = 1, h = 1/2 (dashed line). Bahcall 
& Cen (1992) estimated the cluster mass function using optical data (richness, 
velocities, luminosity function of galaxies in clusters) as well as X-ray data (tem- 
perature distribution function of clusters). Groups poorer than Abell clusters 
have also been included thus extending the mass function to lower masses than 
the richer Abell clusters. Girardi et al. (1998) data are obtained from a sample 
of 152 nearby (z < 0.15) Abell- ACO clusters. As shown, the CDM model that 
does not take account of the non-radial motions over-produces the clusters abun- 
dance. The introduction of non-radial motions (solid line) reduces remarkably 
the abundance of clusters with the result that the model predictions are in good 
agreement with the observational data. This result confirms what found in Del 
Popolo & Gambera (1999) showing how a mass dependent threshold, 5 C (M), 
(dependence caused by the developing of non-radial motions) can solve several 
of SCDM discrepancies with observations. 

The VDF is defined in a similar way to the mass function, namely it is 
the number of objects per unit volume with velocity dispersion larger than cr v . 
Since the velocity dispersion <r v can be observed directly, VDF provides a good 
test of theoretical models. Observed a v comes from the measurement of galaxy 
redshift. The VDF can be calculated starting from the mass function: 

, , ,„,dM 
n(a v ) = n(M)— (6) 

The cumulative VDF can be obtained integrating Eq. @: 

POO 

n(> cr v ) = / n(<7 v /)d<v (7) 



In order to use Eq. (||) to calculate the VDF we need a relation between the 
velocity dispersion, a Y , and mass, M. To determine the relation between <t v and 
M we use both the relation for the typical virial temperature and the result by 
Thomas & Couchman (1992) and that by Evrard (1989, 1990, 1997) found in 
N-body simulations. We find that the necessary relation between a v and M is 
given by: 

(Evrard 1989; Lilje 1990). In Fig. 2 we compare the VDF obtained from a CDM 
model taking account of non-radial motions (solid line) with the CfA data by 
Zabludoff et al. (1993) (full dots) based on their survey of R > 1 Abell clusters 
within z < 0.05 and with the data by Mazure et al. (1996) (full triangles) and 
Fadda et al. (1996) (open squares). Mazure et al. (1996) data are obtained 
from a volume- limited sample of 128 Raco > 1 clusters while that of Fadda et 
al. (1996) are obtained from a sample of 172 nearby galaxy clusters (z < 0.15). 
We also plot the VDF obtained from a CDM model without non-radial motions 
(dashed line). The SCDM model predicts more clusters than the CfA observation 
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Figure 2. Cumulative VDF calculated using a CDM model without 
taking into account non-radial motions (dashed line) and taking into 
account non-radial motions (solid line) compared with Zabludoff et al. 
(1993) data (full dots) and with those by Mazure et al.(1996) (full 
triangles) for R > 1 clusters and Fadda et al. (1996) (open squares) 
for R > —1 clusters. The theoretical curves are obtained using a cr v -M 
relation with zero scatter. 

except at cr v ~ llOOkm/s. As reported by Jing & Fang (1994) the SCDM model 
can be rejected at a very high confidence level (> 6a). The reduction of the 
normalization reduces the formation of clusters, thus resolving the problem of 
too many clusters, but leads to a deficit at a v ~ llOOkm/s. When non-radial 
motions are taken into account (solid line) we obtain a good agreement between 
theoretical predictions and observations. Both CDM and CDM with non-radial 
motions predict more clusters of low velocity dispersion (a v < 300km/s) than 
the observation. This discrepancy is not significant because the data at a v < 
300km/s could be seriuosly underestimated (Zabludoff et al. 1993; Fadda et al. 
1996; Mazure et al. 1996). 



3. Non-radial motions and the shape of clusters 

Most clusters, like elliptical galaxies, are not spherical and their shape is not due 
to rotation (Rood et al. 1972; Dressier 1981). The perturbations that gave rise 
to the formation of clusters of galaxies are alike to have been initially aspherical 
(Barrow & Silk 1981; Bardeen et al. 1986) and asphericities are then amplified 
during gravitational collapse (Icke 1973; Barrow & Silk 1981). The elongations 
are probably due to a velocity anisotropy of the galaxies (Aarseth & Binney 
1978) and according to Binney & Silk (1979) and to Salvador-Sole & Solanes 
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(1993) the elongation of clusters originates in the tidal distortion by neighboring 
protoclusters. In particular Salvador-Sole & Solanes (1993) found that the main 
distorsion on a cluster is produced by the nearest neighboring cluster having 
more than 45 galaxies and the same model can explain the alignement between 
neighboring clusters (Oort 1983; Plionis 1993) and that between clusters and 
their first ranked galaxy (Rhee & Katgert 1987; van Kampen & Rhee 1990; 
West 1994). 

The observational information on the distribution of clusters shapes is sometimes 
conflicting. Rhee et al. (1989) found that most clusters are nearly spherical with 
cllipticities distribution having a peak at e ~ 0.15 while Plionis et al. (1991) 
found that clusters are more elongated with the peak of the ellipicities distribu- 
tion at e ~ 0.5. 

To study the effect of non-radial motions on the shape of clusters we shall use a 
model introduced by Binney & Silk (1979). In that paper they showed that tidal 
interactions between protoclusters and the neighbouring protostructures should 
yield prolate shapes (before virialization) with an axial ratio of protostructures 
of ~ 0.5, the typical value found in clusters. After virialization the pre-existing 
elongation is damped and the axial ratio leads to values of about 0.7-=-0.8, that 
are higher with respect to observations. As observed by Salvador-Sole & Solanes 
(1993) this last discrepancy can be removed taking into account that tidal in- 
teraction keeps going on after virialization and that on average the damping of 
elongations due to violent relaxation is eliminated by its growth after virializa- 
tion. Then this growth restores a value of e near the one that clusters had before 
virialization. 

According to the quoted Binney & Silk (1979) model, an initially spherical 
protostructure (e.g. a protocluster) of mass M having at distance r(t) from its 
centre a series of similar protostructure of mass M' shall be distorted and the 
final ellipticity can be calculated using the following equation: 

e^^y/ 2 G(2, Xl ) (9) 

where G(2, x\) is defined in the quoted paper (see Eq. 13c) and (a. 2 ) 1 / 2 for a 
CDM without non-radial motions is given by (Binney & Silk 1979) while for a 
CDM with non-radial motions we use the following relation: 



(a* ) 



2\l/2 _ 



9^ 



N(M',r)dM'] 1/2 



M + AP 



In Fig. 3 we show as e declines with mass in agreement with the observations. 
Besides, we show as in a CDM model that takes into account non-radial motions 
(solid line) e is smaller than in the simple CDM (dashed line). This show as 
non-radial motions reduce the elongation of clusters. For a cluster of 1O 15 M0 we 
get a value of e ~ 0.5, if non-radial motions are excluded, while e ~ 0.43 when 
non-radial motions are taken into account. Increasing the mass, as expected, 
clusters tend to become more and more spherical. 
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Figure 3. Ellipticity, e, of clusters versus mass, M. The dashed and 
solid lines represent e for a CDM without and with non-radial motions, 
respectively. 

4. Conclusions 

Here, we have studied how non-radial motions change the mass function, the 
VDF and the shape of clusters of galaxies, using the model introduced by Del 
Popolo & Gambera (1998a; 1999). We compared the theoretical mass function 
obtained from the CDM model taking into account non-radial motions with the 
experimental data by Bahcall Sz Cen (1992) and Girardi et al. (1998). The 
VDF was compared with the CfA data by Zabloudoff et al. (1993), by Mazure 
et al. (1996) and by Fadda et al. (1996). Taking account of non-radial motions 
we obtained a notheworthy reduction of the discrepancies between the CDM 
predicted mass function, the VDF and the observations. Non-radial motions are 
also able to change the shape of clusters of galaxies reducing their elongations 
with respect to the prediction of CDM. 
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